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Abstract 

We elaborate on boundary conditions for Ginzburg-Landau (GL) theory in the case of external currents. We 
implement a self-consistent theory within the finite element method (FEM) and present numerical results for a 
two-dimensional rectangular geometry. We emphasize that our approach can in principle also be used for general 
geometries in three-dimensional superconductors. 
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1. Introduction 

There is an increasing interest of using supercon- 
ductors for example in the design of products where 
a reduced size and weight is crucial, such as gen- 
erators in the nacelle of large off shore wind tur- 
bines and electrical motors in azimuth thrusters for 
ships [1]. Modeling of superconductors for techno- 
logical applications is a truly multiscale problem. 
The remarkable development of computers for sci- 
entific calculations makes it practically possible to 
take further steps in the sophistication of quanti- 
tative modeling of superconductors with transport 
currents. FEM modeling of superconducting tapes 
in realistic technological devices can today be per- 
formed under the assumption of an empirical ex- 
pressions for the resistivity [2]. A magnetic field is 
generally expelled from the interior of a supercon- 
ductor. However, in type-II superconductors, quan- 
tized magnetic flux can penetrate from the surface 
and form vortices which increase the macroscopic 
resistivity. The Ginzburg-Landau (GL) theory is a 
celebrated tool for theoretical modelling of super- 
conductors [3-5] such as for example of vortex dy- 



namics in type-II superconductors, see e.g. [6]. The 
focus in the present paper is to discuss implementa- 
tions of GL theory for the microscopic modelling of 
currents through type-II superconductors. 

2. Theoretical formalism 

The time- dependent GL theory is strictly valid 
for gapless superconductors and for temperatures 
close to Tc, but have been confirmed experimentally 
for a wide range of temperatures [7]. We do not 
study the temperature dependence here, moreover 
we fix the electromagnet gauge to be the (incom- 
plete) Weyl gauge, i.e. the electric scalar potential 
is zero cp = througout the paper. The dimension- 
less time-dependent Ginzburg-Landau equations to 
be solved self-consistently are 



atV = -(^V + A)'^ + ^-|^|2v, 
adtA = ^ (^*V^ - ipvr) - 1^1' A 
-V X (V X A - B(°)) , ren; 




Preprint submitted to Physica C 



8 November 2011 



Here, ^p is the order-parameter for the superconduct- 
ing Cooper-pair condensate, while A is the mag- 
netic vector potential. Moreover, B^^^ (r) is the di- 
mensionless applied magnetic field, while the inte- 
gral operator in the lower equation is the current- 
induced magnetic field following Biot-Savarts law, 
where J = J^^^ -h J^^'^ is the sum of super- currents 
and normal- currents [7] 

j(^) = (i/;V7/^* - 7/^*V7/^)-|7/^|^ A; J(^^ = -adtA. 

(2) 

Note that the last term in the lower equation in (1) 
differs from the standard time-dependent GL equa- 
tions, but is instrumental in incorporating trans- 
port currents as we will detail below. The parame- 
ters in the problem are a, />:, / and B^^\ They are 
related to the (unsealed) physical quantities (P) as 
follows: ap = a/ (/io^^^), where a is the conduc- 
tivity of the normal current, /xq is the permeability 
of vacuum, is a phenomenological diffusion coef- 
ficient [5] ; = A/^ is the Ginzburg- Landau parame- 
ter (approximately independent of temperature), it 
is the ratio between the London penetration depth X 
for external magnetic fields and the coherence length 
^ of the (Cooper-pair) condensate. In the present 
formulation we use A as lengthscale, while time is 
scaled according to ^'^/D. The transport current to 
be modeled is Ip = hi/ (/io^<fA), where h is Plancks 
constant divided by 27r and g = 2e is the electric 
charge of a Cooper pair. The magnetic vector po- 
tential is scaled according to Ap = hA/ (q^), hence 
the magnetic field isBp = VxAp = h'B / (g^ A) . 

We now discuss the qualitative behaviour of cur- 
rents within the superconductor in general terms. 
The super-currents J^*^ is zero at the points in con- 
tact with the normal conductors. Super-currents 
carry the main part of the transport current within 
the superconductor, except in the vicinity of vor- 
tices where super-currents encircle the centers of 
the vortices (see Fig. 1). The normal- currents J^^^ 
carry the total current close to the normal conduc- 
tors. Normal-currents are also present at positions 
where the Cooper-pair density is penetrated by 
moving vortices, the normal-currents here are re- 
sponsible for dissipation inside the superconductor 
and build up a macroscopic resistance. Hence from 
a technological point of view, it is important to op- 
timize geometries, dimensions, material properties 
and pinning of vortices, in order to obtain desirable 
current- volt age relations for superconductors. How- 
ever, in general there are rather few articles treating 
GL theory for transport currents, see [8-11] for 



examples. 

Crucial for the modeling of transport is how to set 
up the boundary conditions (BC). From a physical 
point of view, it involves a condition for the super- 
currents J^^^ • n = not to penetrate any boundary, 
a specification of the magnetic field B at the bound- 
ary, and conditions for the normal-currents J^"^^ at 
the boundaries. The vacuum- superconductor BC for 
the standard time-dependent GL equations are [11] 

(^^V + A^ V^.n = 0; VxA = B^^); -cr^tA-n = 0, 

(3) 

while for regions where the current is injected, the 
metal- superconductor BC are 

'0 = 0; V X A = B^^^; -adtA • n = J^^^ • n. (4) 

Here, n is the outgoing normal vector to dfl and 
J^^^ is the external current density. The exter- 
nal magnetic field is the sum of the applied and 
current-induced fields B^^^ = B*^^^ + B*^^^ as shown 
in the lower equation (1). Note that the applied 
field may result from other superconductors, for 
example in a generator or electric motor. In a 
general three-dimensional geometry the implemen- 
tation of the presented theory may be challenging 
and computationally very demanding. Therefore 
we here give results for illustrative purposes of a 
standard geometry in numerical GL simulations, 
the two-dimensional {z = 0) rectangle. That is a 
superconductor of length L^ and width Ly, with 
the center at origo, and in contact with (infinitely 
long) homogeneous normal conductors of width Ly 
at the two boundaries where x = :hLx/2. 

One approximation to the problem corresponds 
to using the standard GL equations and a linear 
function for the magnetic field at the boundary [8] 

B^^) =B^^)-hP^^^^^e^. (5) 

Ly 

This corresponds to using the same current-induced 
magnetic field as that from a superconductor with 
a constant current density. This assumption does 
not take into account the interplay between the dy- 
namics of the currents and the magnetic field at the 
boundary of the superconductor. The true currents 
are larger close to boundaries and vary substantially 
due to the local response of applied fields and corre- 
lates with the vortex dynamics present. In order to 
focus on the role of vortex dynamics, we set B^«) = 
in the following. Let us show here how Eq. (5) is then 
related to the standard BC of Eqs. (3) and (4) and 
that they are redundant. The third BC in Eq. (4) can 
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be reformulated to A = — J^^^ /a • t where the cur- 
rent goes through the boundaries. However, this BC 
needs in fact not to be exphcitly implemented, since 
the following relation between the second and third 
BC in Eq. (4) holds for the left and right boundary 

A = ^-—t ^ V X A = -^y e^. (6) 

(7 Ly Ly 

Indeed, the right form of the above BC gives a cur- 
rent density at the left- and right-boundaries that 
are parallell to e^^, since according to Ampere law 
for the current density (assuming a static electric 
displacement field = 0) 

J{x = ±L,/2, = V X B = ^ e„ (7) 

Ly 

where 

/ = / J^'^-e,c^^, j(-) = ^e,. (8) 

J-Ly/2 

Hence, one can implicitly reproduce the third BC 
of Eq. (4) by a proper choice of the magnetic field 
at the boundaries. Whatever choice is made from 
Eq. (6) for the left and right BC of the rectangle, 
the upper and lower BC should now fulfill 

B^^) (x, y = ±Ly/2) = B^^) = ±1/2. (9) 

We have confirmed in numerical calculations that we 
can equivalently use either of the two possibilities in 
Eq. (6) for the left and right BC, and Eq. (9) for the 
upper and lower BC. 

In order to treat the BC self-consistently, the mag- 
netic field at the boundaries should be calculated 
from the true current density and coupled into the 
GL equations as formulated in Eqs. (1). In the fol- 
lowing we discuss this general problem. We use the 
left part of Eq. (6), i.e. we assume that the non- 
trivial current dynamics inside the superconductor 
do not affect the normal conductors that it is con- 
nected to, and t/j = for the left and right BC. Then 
we solve the lower equation of (1), with A • n = 
and S/tp • n = for the upper and lower boundaries. 
The dynamics of EM-fields in the surrounding half- 
spaces are not considered. 

3. Numerical results 

We use a FEM implementation with cubic La- 
grange shape functions and a time-dependent adap- 
tive mesh refinement [12] to solve Eqs. (1). The to- 
tal transport current through the superconductor 
is monitored numerically with a projection onto Gx 




Fig. 1. Super-currents J^^^ (arrows) of Eq. (2) for; (a) a 
self-consistent-; and (b) a linear-; calculation of the curren- 
t-induced magnetic field. Background gray-scale shows the 
corresponding density \tp\'^ of Cooper-pairs. Parameters used 
8,re a = K = 4 (type-II), / = 1.55 and B(^) = 0. The size 
of the rectangle is Lx = Ly = lOA. The dimensionless time 
shown here is t = 200, however, the external current / and 
the amplitude of the self-consistent BC are applied gradually 
to keep it numerically tractable. At t = 200 it is 8% of the 
self-consistent amplitude that is mixed into the BC, which 
is already sufficient to demonstrate qualitative differences. 
Convergence has been confirmed with different meshes. 

(direction of transport) with integration along lines 
orthogonal to (i.e. I = J Jx{x^y)dy). 

We give examples of snapshots of the super- 
currents and vortex structure for a self-consistent 
implementation of Eqs. (1) in Fig. 1 (a), and for a 
simulation, using the approximate BC of Eq. (5), 
in Fig. 1 (b). On the boundary r G we use a 
simplified notation for the magnetic field at the 
upper- /lower- (U/L) boundaries of the rectangle 

(10) 

In practice we implement a numerical integration 
over the superconductor for each time-step. For the 
trivial rectangular geometry, for example 

(11) 
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Fig. 2. Current-induced magnetic fields at the upper and 
lower boundaries corresponding to Fig. 1 (a), solid curves 
shows results of Eqs. (10)-(12). Dashed lines shows the ap- 
proximate theory of Eq. (9). Parameters as in Fig. 1. 

is the contribution to Bu from the superconductor. 
Then B^^^^ is obtained by an analogue construc- 
tion. The contribution B^^^ from the integration 
over the two semi-infinite normal conductors, where 
j(^) = I /Ly is an even function in odd in 



(NC) 



such that B^ 
lated only once 



B^^^ , and needs to be calcu- 



(NC) 



+atanhl 



+atanh| 



(12) 



In Fig. 2 Bjj/L from Eqs. (10)-(12) is compared 
with the linear approximate theory of Eq. (5). In 
the numerical time evolution we observe the cre- 
ation of antivortex- vortex pairs symmetrically at the 
upper and lower boundaries (Fig. 1), followed by 
annihilation within the superconductor when they 
meet at the horizontal axis. Note in particular that 
Bl = —Bu here (Fig. 2) due to the symmetry of 
antivortex-vortex pairs in the absence of an exter- 
nal applied magnetic field. If there is also an applied 
field B*^^^ 7^ 0, vortices are created at the boundary 
with highest external magnetic field (if above a crit- 
ical value). Generally then doesn't obey any 
symmetries but refiects the non-trivial dynamics of 
external fields and vortices present. 

Following the positions of vortices and current- 
induced fields through time, one can qualitativ- 
elly explain observed differences between the sim- 
ulations using self-consistent vs linear BC [i.e. 
Figs. 1 (a) vs (b)]. We here refer to those as sim- 
ulations (a) and (b), respectively. For example, 
the three vs two antivortex-vortex pairs present at 
t = 200 in Figs. 1 (a) and (b), can be explained as 
a memory effect for t ~ 145 (not shown) when vor- 
tices communicates with the BC via Eq. (11). At 
this point in time, both simulations are similar, due 
to the gradual increase of the amplitude for the self- 
consistent BC, and show three antivortex-vortex 



pairs. For the (b) simulation, two more pairs then 
enter such that a short time before the three central 
pairs annihilate, there are five pairs present simul- 
taneously. However, for simulation (a), the three 
pairs present at t ~ 145 in effect reduce the current- 
induced fields, since a vortex introduces an opposite 
directed (to e^^) super-current with a smaller dis- 
tance to the (e.g.) y = Ly/2 boundary compared to 
the lower part of the same vortex, and hence locally 
lower the value of Eq. (11). New antivortex-vortex 
pairs in simulation (a) are therefore delayed in time 
until older pairs have been annihilated. 

4. Summary 

We have presented self-consistent boundary- 
conditions to model transport currents with the 
time-dependent Ginzburg-Landau equation. Nu- 
merical results for the vortex dynamics confirm 
qualitative differences compared to the standard 
GL equation and BC, however, a practical drawback 
is a substantially higher computational cost. 
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